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Abstract: The complete mitochondrial DNA D-loop sequences from 10 stud Brahman cattle were sequenced and ana- 
lyzed. The results showed that the genetic diversity of Brahman cattle was rich; the rate of nucleotide variation, haplotype 
diversity and nucleotide diversity were 6.25% , 0.978 + 0.054 and 0.014 30 + 0.008 68, respectively. Nine haplotypes 
were defined and fell into two distinct lineages, suggesting that Brahman cattle have both Bos indicus (Zebu) and B. tau- 
rus genetic background. The taurine haplotypes were predominant at 90% and only Brah-6 belonged to the Asian zebu mt- 
haplotype. This indicates that Brahman cattle was one of the zebu breeds and inherited the excellent characteristics of both 
the Asian zebu and European beef cattle, such as easy calf delivery, high quality beef, heat tolerance and resistance to 
various parasites. Breeders introduced Brahman cattle to improve the productivity and adaptability of native cattle. The Ze- 
bu has evidently frequently introgressed into the modern taurine breeds. As for modern zebu breeds, B. taurus also highly 
contributed to their formation, except for the Asian zebu. Furthermore our results also confirm the hypothesis that B. indi- 
cus has undergone a separate domestication event and originated from the Indian subcontinent. 


Key words: Brahman cattle; mtDNA D-loop polymorphism; genetic background 


婆罗 门 牛 mtDNA D-loop 变异 及 其 遗传 背景 
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(1. 云南 省 肉牛 和 牧草 研究 中 心 , 云南 昆明 650212; 

2. 中 国 科学 院 昆 明 动 物 研 究 所 云南 省 畜 禽 分 子 生物 学 重点 实验 室 、 分 子 进 化 与 基因 组 多 样 性 实验 室 , 云南 昆明 650223) 


摘要 : 对 10 KRAVE IA mtDNA D-loop 全 序列 912 bp 测序 , 婆罗 门 牛 遗传 多 样 性 丰富 , 检测 到 的 9 种 单 
倍 型 兼 有 瘤 牛 (了 indicus ) 与 普通 牛 (B. taurus ) 的 遗传 背景 ， 核 苷 酸 变 异 率 为 6. 25% , 单 倍 型 多 态度 为 0.978 + 
0.054, HRA AEA 0.014 30 + 0.008 68。 所 有 单 倍 型 聚 为 明显 的 两 大 分 支 , 婆罗 门 牛 的 大 部 分 单 倍 型 为 普通 
牛 单 倍 型 类 群 , 并 占 绝对 优势 (900% ), 仅 Brah-6 与 亚洲 瘤 牛 聚 在 一 起 , 属于 亚洲 瘤 牛 线粒体 单 倍 型 , 表明 婆罗 门 
TRAE REAT MAET ERRET ATE T AREF, 耐 热 与 体 表 寄生 虫 等 ) 的 瘤 牛 品种 之 

育种 学 家 引种 瘤 牛 的 目的 是 改善 当地 牛 的 生产 力 与 适应 性 , 现代 普通 牛 表 现 出 明显 又 普遍 的 瘤 牛 渐 渗 现象 。 
对 现代 的 瘙 牛 品种 而 言 ， 除 亚洲 瘙 牛 品种 外 ,普通 后 对 其 他 瘤 牛 品种 育成 的 贡献 同样 高 。 KERE (B. indias) 
为 独立 驯化 、 起 源 于 印度 次 大 陆 的 假说 。 
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Domesticated cattle have played a key role in agri- 
culture, economy, culture and religion during the his- 
toric and modern societies of human civilization. Zoolo- 
gists usually classified cattle as two separate species: 
B. indicus, indicine, humped; and B. taurus, tau- 
rine, humpless. However, the complete inter-fertility 
between B. indicus and B. taurus has given rise to 
their consideration as subspecies (Grigson, 1980; Ep- 
stein & Mason, 1984; Payne, 1991), and that all cat- 
tle breeds derived from the now extinct wild ox or au- 
rochs, B. primigenius. It is well known that the first 
domestication of cattle was some 8 000 — 10 000 years 
ago in the Near East. Cattle subsequently dispersed 
around the world and gave rise to all the modern cattle 
breeds, excluding Bali and Mithan cattle ( Grigson, 
1980; Epstein & Mason, 1984; Payne, 1991; Loftus 
et al, 1999). Loftus et al (1994a, b) hypothesized 
that Asian zebu had undergone an independent domesti- 
cation. mtDNA has an evolutionary rate 5 — 10 times 
faster than typical nucleus genes. In addition, the mi- 
tochondrial DNA control region (mtDNA D-loop ) has 
an evolutionary rate which is faster than other mitochon- 
drial regions. MtDNA was used to clarify the origin, 
domestication and genetic background of domesticated 
cattle. 

Brahman cattle originated from India about 4 000 
years ago. As one of the zebu breeds, it is character- 
ized by long floppy ears and a large hump over the top 
of the shoulder and neck. It has highly developed sweat 
glands and a dewlap and has good heat tolerance and 
parasite resistance. They can travel long distances to 
water and therefore have a larger grazing range. Brah- 
man cattle seldom contract eye cancer and pinkeye. 
They have a considerable reputation for rapid growth 
and a high dressing percentage, or the percentage of 
animals killed for meat (Oklahoma State University 
Board of Regents, 2000). Based on the above, Brah- 
man cattle were introduced to improve the adaptability 
and production of local cattle breeds in the mid-18th 
century. Later, breeders crossbred some beef cattle 
breeds, including the Beefmaster (1/2 Brahman x 1/4 
Shorthorn x 1/4 Hereford) , Santa Gertrudis (3/8 Brah- 
man x 5/8 Shorthorn), Braford (3/8 Brahman x 5/8 
Hereford). In Yunnan province, China, Brahman cat- 
tle were also introduced to cross with Yunnan native 
yellow cattle in the early 1980s. The aim was to cross- 
breed a new beef cattle breed (or strains) that could 
adapt well to the heat and humid weather in South Chi- 
na (Wen et al, 2003). The resulting breed (1/2 Brah- 


man x 1/4 Murray Gray x 1/4 Yunnan Yellow cattle, 
BMY) is breeding inter se. To find out the genetic di- 
versity and background of Brahman cattle, we se- 
quenced the complete mtDNA D-loop. This will cast 
light on more genetic information about the conservation 
of breeds in the tropical and subtropical areas of South 


China. 
1 Materials and Methods 
1.1 Materials 


The blood from ten Brahman cattle was collected 
from the tail artery at Xiaoshao Farm of the Yunnan 
Beef Cattle and Pasture Research Center. 

1.2 DNA extraction and sequencing 

Total DNA was isolated from fresh blood samples 
as described previously (Yuan et al, 2001). D-loops 
were amplified on a Perlin-Elmer DNA thermocycler. 
Primers were designed from the published bovine com- 
plete mtDNA sequence ( V00654, Anderson et al, 
1982). The forward primer was L15718: 5'-CTAA- 
GACTCAAGGAAGA AACTGC-3' and the reverse primer 
was H415: 5'-GACTCATCTAGGCATTTTCA-3'. The 
PCRs were carried out using a standard PCR program 
with 5 min denaturizing at 95 °C , 35 cycles for 1 min at 
94 °C, 1 min annealing at 58 °C, 1 min extension at 72 
°C., and final extension for 10 min at 72 ^C. To ensure 
the accuracy of the sequences, we also used another 
two primers when sequencing reactions were done, 
H16320: 5’-TTAAGAGGAAAGAATGGAC-3’ and 
L16230: 5'-ACCGTGGGGGTCGCTAT-3’. Purification 
of PCR products and sequencing reactions followed Qu 
et al (2003) . 

1.3 Data analysis 

Data analysis was processed according to Qu et al 
(2003). The NJ tree was constructed using the software 
Mega 2.1 (Kumar et al, 1993) . 


2 Results and Analyses 


2.1 Variation in the mtDNA D-loop of Brahman 

cattle 

Nine hundred and twelve bp complete mitochon- 
drial DNA D-loops of ten Brahman cattle was se- 
quenced and 57 polymorphic sites were found (Gen- 
Bank accession numbers DQ887760 — DQ887768 ) . 
There was a variation in length of a 6 — 8 bp poly (C) 
tract found at 764 bp site. The variance rate of nu- 
cleotides was 6.25%. Three of the 57 substitutions 
were due to transversions, reflecting a heavy transition- 


al bias, which was congruent with previous bovine con- 
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trol region variations (Loftus et al, 1994a; Bradley et 
al, 1996; Mannen et al, 1998). The three transver- 
sions were at 150 bp, 162 bp and 831 bp. There were 
also four 1 — bp insertion/deletions at 352 bp, 410 bp, 
770 bp and 771 bp, respectively (Fig. 1). 
cattle showed high genetic diversity, of which haplotype 


Brahman 


diversity, nucleotide diversity and nucleotide difference 
were 0.978 + 0.054, 0.014 30 + 0.00 868 and 13.000 
+ 8.232, respectively. The Tajima's D value was — 
1.4982 (P >0.10), which was in accordance with the 
neutral mutation hypothesis. By aligning the sequences 
nine haplotypes were defined, of which Brah-6 be- 
longed to the classical Asian zebu mt-haplotype and the 
remaining eight haplotypes clustered into the B. taurus 
clade (Bradley et al, 1996; Lai et al, 2005). Brah-1 
occurred in duplicate, and the others occurred once 
each. 

The accession number and cattle breeds of partial 
haplotypes from GenBank and 10 Brahman cattle are 


listed. Jl represents for Japanese black cattle, 


U87633; QC-06-A1 for Qinchuan cattle, AF514784; 


V00654 for Holstein cow; FR482, FR521, FR479 for 
Friesian, AY378146, AY378144, AY378145; CH-04- 
B1 for Chinese Holstein, AF516713; Lulu-1, Lulu-3, 
Lulu-4 for Lulu cattle (Nepal ), AB085918, 
AB085920, AB085921; Hariana-1, Hariana-2 for Har- 
iana (India), AB085922, AB085923; ON2770, 
ON2771，ON2774 for Ongole (India), AY378134, 
AY378135, AY378136; YY-03-G1 for Yueyang Yel- 
low cattle, AY119678; Nellore/Nelore for Nellore 
(Brazil), NC- 005971; TH982 for Tharparkar (Indi- 
a), AY318137; Zwergzebu for German small-typed ze- 
bu, AF492350; Brah-1 - 9 for Brahman sequenced in 
this study, respectively . 
2.2 The phylogenetic status of Brahman cattle be- 
tween B. indicus and B. taurus 

Together with a range of complete D-loop se- 
quences of the genus Bos published in GenBank, we 
aligned the polymorphic sites (Fig. 1). Al the haplo- 
types were clustered into two distinct clades , the 
B. taurus and B. indicus clades (Fig. 2). The 


Brahman cattle carried both B. taurus and B. indicus 


1111122 2222222222 2223333333 3333333333 3344444444 4455555556 7777777777 78888 
2253566602 3455566678 9990011222 2333344555 7901344556 6701112585 1123577889 90345 
6755002838 1918956763 1342818256 8014967026 4450912780 5440126674 4718401238 94160 


V00654 TATTATCGAT GATCCAGCAT GCTGTGTTTG TGTATTTTAT TGG-ACCCCA TGTACGTGCT AAACT—-TTC TCATG 
FR482 人 Ku a dt eh ip eene T ns a 
Brah-5' Ow une i de shaven Gun uiu 和 Uses Ae ens 
Brahe ilv iade EE C ete Ex apu Ae Seed Gs ed CO aa his 
Brah-4 > coiled gad)’ a Sep easiest Ge aide tase WO INA AL rs mede us A l.ss eos [c E 
Brahe?! -oreli aE abe RENE eta a ET ean dex | ecc pL dis ees iG? ket. 
X7503-63- Li ahi Gy wr keen rd GC PET jT a a 人 
LU so 人 Cu Liu AL QUT d te aad ee ee Geeta va damen 
Braáh-l. a a HIR lOS EP DOS C.. zh AUT Grias DU prea wikis 
Lulu-l — ....... Asa ote eed e aut badd uy ni RUN WesCo escrita sinet Tis 66.99 s nnus 
CH-04-Bl ikiia ni Dry Cd A Ac uk SES a [Ou PEE 
QC-06-A1 ........us ee [oe HL Ao hob eee SP td es Ae ee eis Helis Arego a 人 
本 | Cr AIV eds ake ae Kov ende el A sno (sc S ERRA 
Brah-9 o .......... Ke Pd uerus Cl uni obi Re 人 (vu iS. ceu 
Brahem xmv ED Tei es C cece 人 [LT CeCe. ortu 
下 和 Ao Tt onis E CLS der y ONE De ESAME C...T Gs ee Cac cies 
ERSA, o sexe she bescs ACT s Ts rr decide Aem iss eed esed Mos Gass CAs ce vss 
Brah=8 o- llus Aisa Petar hg ko de e anis j M perder Coi stes da Ea! atas 
Brah-6 G.. TCGAG. A. . T. . AT. C AT... ACCCA CACGCCC. -C .. AAG.. TT. ... GT.. A. C 6.6. CC-CCT C... CA 
Nellore G... CGAG. A.. T. . AT. C AT... ACCCA CAC. CCC.-C .. AAG.. TT. ... GT.. A. C. 6.6. CC-CCT C... CA 
YY-03-61 G...CGAG. A.. T. . AT. C AT... ACCCA CAC. CCC. -C .. AAG. . TT. ...GT.. A. C G.G.C--CCT C... CA 
TH982 .G...CGAG. A.. T. AT. C AT. . CACCCA CAC. CCC. C .. AAG. . TT. ...6T.. A.C 6.6. CC-CCT C..CA 
Hariana-2 .G...CGAG. A.. T. . AT. C AT... ACCCA CAC. CCC.-C .. AAG. . TT. .... T.. A. C 6.G. CCC. CT C... CA 
0N2770 G...CGAG. A.. T. . AT... AT... ACCCA CAC. CCC. -C .. AAG.. TT. ...G6T.. A.C 6.6. CC-CCT C. TCA 
0N2774 G... CGAG. A.. TT. AT... AT... ACCCA CAC. CCCC-C ..AAG.. TT. ...GT.. A.C 6.6, C—CCT C..CA 
Lulu-4 . 6G... CGAG. A... T.. AT. C A. C.. ACCCA CAC. CCCC-C .. AAG. TTT. ....T.. ATC 6. 6. CC-CCT C... CA 
Zwergzebu .G...CGAG. A.. T. . AT. C A. C. . ACCCA CAC. CCCC-C .. AAG. TTT. ...GT.. ATC 6... CC-CCT C... CA 
Hariana-1 .6. C. CGAC.. A..... ATGC A. C.. ACCCA CAC. CCCC-C .. AAG. TTT. ...GT.. ATC 6... CC-CCT C... CA 
0N2771 . GC. . CGAG. AG. TT. AT. C A. C. . ACCCA CAC. CCCC-C C. AAG. TTT. ...GT..ATCG .G.CC-CCT CT. CA 
Fig. 1 The polymorphic sites of the mtDNA D-loop of Brahman cattle defined according to 


Anderson et al (1982) and published in GenBank 


wee € o» 


indicates the same base pair, 


indicates deletion/insertion . 
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Fig. 2 A Neighbour Joined tree for Brahman cattle and 
GenBank data, constructed by Mega 2.1 


mt-haplotypes, suggesting that Brahman cattle con- 
tained both indicine and taurine genetic background. 
Ninety percent of the mt-haplotypes were B. taurus . 
The eight mt-haplotypes of Brahman cattle were inter- 
mingled with the Lulu, Qinchun, Xizhen, Holstein 
and Japanese Brown branch. Only Brah-6 mt-haplotype 
belonged to the Asian zebu group, clustering with the 
Indian subcontinental zebu distributed in Nepal, India 
and Pakistan. However, zebu introgression occurred 
frequently; Lulu, Yueyang cattle and other taurine cat- 
tle breeds, showing typical morphological features of 


B. taurus, possessed B. indicus mt-haplotypes . 


3 Discussion 


The Brahman cattle showed high genetic diversity , 
of which haplotype diversity and nucleotide diversity 
were 0.978 and 1.43% respectively, congruent with 
previous results (Expected heterogeneity, Hep = 
0.7396, Qu et al, 2006). Brahman cattle carried both 
B. indicus and B. 


which is one of the leading reasons that Brahman cattle 


taurus mitochondrial lineage, 


possess high genetic diversity (Loftus et al, 1999). 
Kikkawa et al (1995) and Meirelles et al (1999) found 
that nearly all Brahman cows that had no B. indicus 
mt-type resulted from male-mediated introgression. 
This may indicate that Brahman cattle is not a pure ze- 
bu breed because its taurine mt-haplotypes were 90% 
of the total. After the 1960’ s introduction of Indian 
Ongole cattle to Brazilian Nellore cattle, Nellore cattle 


were bred and as expected possessed Asian B. indicus 


mt-type, but B. taurus mt-types were still predominant 
at 79% (Meirelles et al, 1999; Oklahoma State Uni- 
versity Board of Regents, 2000). Although the Butana 
and Kenana cattle in East Africa are characterized by 
zebu features, there was no evidence of the dispersal of 
African zebu with its hybrid origin in the Asian zebu 
clade (Loftus et al, 1994b). Tt seems that many peo- 
ple have underestimated the contribution of B. taurus 
during the formation of the zebu breeds. Modern zebu 
breeds have both B. taurus and B. indicus mt-types 
and they were bred through male-mediating and grad- 
ing-crossing. The newly bred beef cattle breeds were 
improved in adaptability and productivity compared 
with the local breeds (Meirelles et al, 1999). Thus 
Brahman cattle were originally selected from the Indian 
subcontinent and spread widely in the tropical and sub- 
tropical areas, or hot and dry environments, of Ameri- 
ca and Africa. Mongolian cattle, with B. taurus mor- 
phological features, carried distinct B. indicus mito- 
chondrial lineages (2096) and they didn't inherit the 
selective advantages of zebu cattle, such as heat toler- 
ance, tick resistance and a distinct hump, but are well 
adapted to the northern cooler weather as northern cat- 
tle (Mannen et al, 2004). 


have hybrid origins and have undergone zebu introgres- 


Most cattle are known to 


sion, such as Philippine native cattle (Kikkawa et al, 
1995), Lulu (Fujise et al, 2003) , Yunnan cattle (Yu 
et al, 1999), Qinchuan, Jinnan (Lei et al, 2004), 
Guanling (Liu et al, 2005) and Sichuan cattle (Lai et 
al, 2005). Through crossbreeding the BMY cattle con- 
tained both taurine and indicine matrilineages (Data 
not shown) just as Yunnan cattle did. The BMY breed 
is well adapted to the tropical and humid weather in 
South China and can still graze in snow. These eco- 
nomic characteristics have been stably inherited up to 
now (Wen et al, 2003; Yunnan Beef Cattle and Pas- 
ture Research Center, 2005 - 2006) . 

Bos indicus and B. taurus were clustered into two 
distinct clades (Fig. 2). Brah-6 was the closest to the 
Asian zebu and the other eight haplotypes were inter- 
mingled into the B. taurus clade, including 
Qinchuan, Xizhen, Holstein and Japanese Brown (Lof- 
tus et al, 1994a, b; Bradley et al, 1996; Mannen et 
al, 1998; Troy et al, 2001; Magee et al, 2002; Fu- 
jise et al, 2003; Lai et al, 2005). It is worth noting 
that there is a predominant taurine genetic background 
in most of today’s zebu breeds except for the Asian ze- 


bu; B. taurus contributed more to today’s zebu breeds 
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than the Asian zebu did to B. taurus. All the B. indi- 
cus mt-haplotypes belonged to the Asian zebu lineage, 
and the Indian Hariana-2 was the proximal breed of the 
B. indicus clade, indicating that Hariana-2 may be 
from one of the original breeds. We also found strong 
support that the zebu has undergone an independent 
domestication and probably originated from the Indian 
subcontinent (Loftus et al, 1994a; Bradley et al, 
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